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Abstract. - Starting with a scalar field in a thermal bath and using the one loop quantum 
correction potential, we rewrite the Klein-Gordon equation in its thermodynamical representation 
and study the behavior of this scalar field due to temperature variations in the equations of motion. 
We find the generalization of a Gross-Pitaevskii like equation for a relativistic Bose gas with finite 
temperature, the corresponding thermodynamic and viscosity expressions, and an expression for 
the postulate of the first law of the thermodynamics for this BECs. We also propose that the 
equations obtained might help to explain at some level the phase transition of a Bose-Einstein 
Condensate in terms of quantum field theory in a simple way. 



Phase transitions are changes of state, related with 

■ changes of symmetries of the system. The analisys of 

■ symmetry breaking (SB) mechanisms have turn out to 
, be very helpful in the study of phenomena associated to 
' phase transitions in almost all areas of physics. Bose- 
Einstein Condensation (BEC) is one topic of interest that 
uses in an extensive way the SB mechanisms [3 , its phase 
transition associated with the condensation of atoms 
in the state of lowest energy and is the consequence of 
quantum, statistical and thermodynamical effects. The 

■ concept of BEC has in recent years emerged in an array 
of exciting new experimental and theoretical systems and 
is now a common phenomenon occurring in physics at 
all scales, from condensed matter to nuclear, elementary 
particle physics, astrophysics and now cosmology, with 
ideas flowing across boundaries between fields. The net 
results of these varied interests in this phenomenon has 
done that the emphasis gradually shifted from the study 
of nonrelativistic Bose systems to that of relativistic ones, 
[3] . The increasing studies in the correspondence between 
SB in condensed matter physics and cosmology are one 
of the motivations of this research |SHZ]. Although phase 
transitions happen in condensed matter systems and the 
Universe, their status is not similar in both cases. For 
condensed matter, more transitions are known, for which 
our knowledge is better [SHTO]. 



ory raise important challenges about their possible man- 
ifestation in condensed matter systems, which can be 
drawn by non-relativistic or relativistic scalar field (SF) 
theories in the framework of the Ginzburg-Landau theory 
of condensed systems. By investigating the massive Klein- 
Gordon equation we can be able, in principle, to simulate 
a condensed matter system. Since many particles have 
mass, this can be an essential step in building realistic 
analogue models. Here, our aim is to study a relativistic 
model made up of a real SF together with the possibil- 
ity that this SF might undergo a phase transition as the 
temperature of the system is lowered. To analyze the SB 
undergone by the system we take a model having an U(l) 
symmetry, considering the easiest case of a double- well in- 
teracting potential (Mexican-hat potential) for a real SF 
$(a;, t) that goes as 
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The results from finite temperature quantum field the- 



Here one important idea, to which we shall refer, is that 
of identifying the order parameter ^' which characterizes 
the phase transition with the value of the real scalar 
quantum field $. 

From quantum field theory we know that the dynamics 
of a SF is governed by the Klein-Gordon (KG) equation, 
it is the equation of motion of a field composed of spin- 
less particles. In this case we will add a first order auto- 
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interaction potential to the SF such that the KG equation 
will be given by 



In terms of function '5, the KG equation ^ reads, 
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where the interaction potential will be denoted by (f>, and 
as will be seen later on it represents the gravitational po- 
tential, as usual the D'Alambertian operator is. 
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Here we will be interested in studying the properties, be- 
havior and interactions of the SF with other particles, all 
of them interacting inside a thermal bath in a reservoir 
that can have interaction with its surroundings. In this 
case, the SF can be described by the potential ([T]) extended 
to one loop, see [TTM13I . In this case the finite-temperature 
effective potential denotes the free-energy density associ- 
ated with the field $, and goes as: 
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where A = \/{fi^c^) is the parameter describing the 
interaction, "fh? = rr?(? jf? is the scalar particle mass, ks 
is Boltzmann's constant, h is Planck's constant, c is the 
speed of light and T is the temperature of the thermal 
bath, this result includes both quantum and thermal 
contributions. 

From ([S]), the energy density can be computed, 
= + Ts*, where = -dVT{^)ldT. When the 

temperature T is high enough, one of the minimums of 
the potential ([3]) is $ = 0. At this point, the SF density 
is equal to ~ ToT^, with Tq = Tr^fcl/SO^i^^^ We 
suppose that the temperature is sufficiently small so that 
the interaction between the SF and the rest of matter 
decouples. After this moment, the term TpT"* is not 
longer important, as for sufficiently low T the term that 
goes as can be dropped out. 

The critical temperature where the minimum of the po- 
tential $ = becomes a maximum and at which the sym- 
metry is broken is. 



Now for the SF we perform the transformation 
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where k is the scale of the system, which is to be de- 
termined by a experiment, from now on we will take c — 1. 
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where we have kept just the equation for the real part, 
the complex conjugate can be described in the same way. 
The notation used is; ' = d/dt and = = p. © 
is KG's equation rewritten in terms of the function ^ 
and temperature T. This equation is an exact equation 
defining the field 4'(x,t), where (p defines the external 
potential acting on the system and the terms in A 
represent the interaction potential within the system. 

From the mathematical point of view, is a type of 
nonlinear Schrodinger equation. From now on equation 
([S]) will be considered as a generalization of the Gross- 
Pitaevskii equation for finite temperatures and relativistic 
particles. This is because at second order, when T = 
and in the non-relativistic limit, — >• V^, eq. ([5]) 
becomes the Schrodinger equation with an extra term 
In this limit this is known as the Gross-Pitaevskii 
equation for Bose- Einstein Condensates (EEC), i.e, an 
approximate equation for the mean-field order parameter 
in the classical theory, [M]. The static limit of equation 
([5]) is known as the Ginzburg-Landau equation. 

In what follows we transform the generalized Gross- 
Pitaevskii equation ([5]) into its analogous hydrodynamical 
version, [15l[l6], for this purpose the ensemble wave func- 
tion ^/ will be represented in terms of a modulus p and a 
phase S as. 
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where the phase S{x,t) is taken as a real function. 
As usual this phase will define the velocity. Here we 
will interpret p{x,t) = p/Mr as the rate between the 
number density of particles in the condensed state, 
p = mriQ = mNg/L^, being A'o the number of particles in 
condensed state and Mt the total mass of the particles 
in the system, being both, S and p, functions of time and 
position. The concept of SB is often used as a sufficient 
condition for BEG. The assumption that the ground 
state can be macroscopically occupied, is nothing but 
Einstein's criterion for condensation. 

So, from this interpretation we have that when the 
KG equation oscillates around the $ = minimum, the 
number of particles in the ground state is zero, p = = p. 
Below the critical temperature Tc, close to the second 
minimum, <i>^„j„ = kg{T^ - T^)/4:, the density wih 



oscillate around p 
equation 



Kfk%{T^ - T2)/4 as can be seen by 



From ([5]) and ([6]), separating real and imaginary parts 
we have 



p + —{pU'S VSVp - pS) = 0, 
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and 
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System Q is completely equivalent to equation ([5]) and 
is the hydrodynamical representation of the latter. 

To simplify system ^ from now on we will neglect sec- 
ond order time derivatives and products of time deriva- 
tives. In this limit we arrive to the non-relativistic system 
of equations (jH]), which reads 



Taking the gradient of ([7]) and using the definition 



V = —VS 
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then the velocity can be determined by the scalar function 
S, with this we obtain, 



p + V-ipv)~-{pS) 
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Notice that in (j9bl) h enters on the right-hand side through 
the term containing the gradient of p. This term is usu- 
ally called the 'quantum pressure' and is a direct conse- 
quence of the Heisenberg uncertainty principle, it reveals 
the importance of quantum effects in interacting gases. 
Multiplying by p, (j9b| can be written as: 



pi) + p{v ■ V)v = pF^ -Vp + pFq + V • (t, 



(10) 



where — ~'V(j> is the force associated to the external 
potential (j), p can be seen as the pressure of the SF gas 
that satisfies the equation of state p = wf?, Vp are forces 
due to the gradients of pressure and uj — —'iX/Am^ is an 
interaction parameter, Fq — —^Uq is the quantum force 
associated to the quantum potential, |17j . 



Uq = - 

and V • (7 is defined as 
V • cr = —pvS - 
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where the coefficient C, is given by 
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and the term V(ln/5) can be written as 



V(lnp) = -V(V • v) - V[V(lnp) • v] + —V[^{pS)] 



p + V- (pv) 
pi) + p{v ■ V)i; 
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Vp + pFq + V -a. 
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Equation p2ap is nothing but the continuity equation, 
and (jl2bp is the equation for the momentum it contains 
forces due to the external potential, to the gradient of 
the pressure, viscous forces due to the interactions of the 
condensate and finally forces due to the quantum nature 
of the equations. 

In this case V(ln/5) now reads 



Thus 



V(lnpj = -V(V • v) - V[V(lnp) 
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V • cr = -—k^pTVT - C [V(V • v) + V[V(ln/3) 
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Here we interpret the function V • a as a viscosity 
expression, it contains terms which are gradients of the 
temperature and of the divergence of the velocity and 
density (dissipative contributions). The measurement of 
the temperature dependence in this thermodynamical 
quantity at the phase transitions might reveal important 
information about the behavior of the gas due to particle 
interaction. 

In what follows we will derive the thermodynamical 
equations from the hydrodynamical representation. We 
can derive a conservation equation for a function a, start- 
ing with the relationship 



{pa) — pa + ap 



(13) 



where a can take the values of and [/g, both of them 
fulfil equation p^ . Using the continuity equation (jl2ap 
in ([T5|) we obtain, 

[pa) -\- V • {pva) = -pv ■ Fa + pa- 

Nevertheless, this procedure is not possible for a because 
in general we do not know it explicitly, only in some cases 
it might be possible to integrate it. 

Observe how the quantum potential Uq also fulfills the 
following relationship 



piiQ + v ■ (pvp) = 0, 



(14) 
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which follows by direct calculation, and where we have 
defined the velocity density Vp by 

(Vlnp); 



Am? 



which can be interpreted as a velocity flux due to the po- 
tential Uq. Using the continuity equation ()12a|) . equation 
(fTl|) can be rewritten as 

{pUq) + V • [pUqv + Jp) + pvFQ=0 (15) 

where we have defined the quantum density flux 

Jp= pVp. 

Equation (jl5p is another expression for the continuity 
equation of the quantum potential Uq. 

As we know, in general (for non-relativistic systems), 
the total energy density of the system e is the sum of the 
kinetic, potential and internal energies |TS], in this case 
we have an extra term Uq due to the quantum potential, 

e = /oe = ipw^ + pcj) + pu + pUq (16) 

and where u is the inner energy of the system. Then from 
(|16p we have that u will satisfy the equation 



{pu)-+V ■ Ju-V ■ Jp 



(17) 



being the energy current, given by a energy flux and a 
heat flux, Jq, 

Ju = PUV + Jq~ pV, 

where 'V ■ Jq = v(V • cr), expression that as we can see is 
related in a direct way to the velocity and gradients of tem- 
perature in the condensate, and is the one that shows in 
an explicit way the temperature dependence of the ther- 
modynamical equations. With these definitions at hand 
we have, 



(pu) ' + V • {pVU + Jq - pV - J p) 



-pV-v. (18) 



In order to find the thermodynamical quantities of the 
system in equilibrium (taking p as constant on a volume 
L), we restrict the system to the regime where the auto- 
interacting potential is constant in time, with this condi- 
tions at hand for ((T5)) we have: 

ipuy+ V • {pvu + Jq-pv- Jp) = -pV ■ V (19) 

From (|19p we can have a straightforward interpretation of 
the terms involved in the phase transition. As always the 
first term will represent the change in the internal energy 
of the system, — pV • v is the work done by the pressure 
and V • V is related to the change in the volume, Jq con- 
tains terms related to the heat generated by gradients of 
the temperature VT and dissipative forces due to viscous 



forces ^ V(V • t>) and finally but most important we have 
an extra term, V • Jp, due to gradients of the quantum 
potential ([TT|). 

Integrating this resulting expression on a close region, 
we obtain 
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J pu dV + <j) {J q + pv) ■ n dS — <j) Jp-ndS 



Equation ()19|) is the continuity equation for the internal 
energy of the system and as usual, from here we have 
an expression that would describe the thermodynamics of 
the system in an analogous way as does the first law of 
thermodynamics, in this case for the KG equation or a 
BEG. This reads 



dU = dQ + dAQ - pdV 



(20) 



where as always, U = ^ pu dV is the internal energy of 
the system, [T9|, and as we can see, its change is due to 
a combination of heat Q added to the system and work 
done on the system (pressure dependent), and 



AAq 
dt 
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is the corresponding quantum heat fiux due to the 
quantum nature of the KG equation. The second term 
on the right hand side of equation ([20| would make the 
crucial difference between a classical and a quantum first 
law of thermodynamics. 

From hereafter we study the transition between the 
$ = state to the low energy one with T < Tc close to 
the minimum $mm = ks^/T^ — 7^/2. 

During the time when T » Tc there are not scalar 
particles in the ground state. We will suppose that the 
scalar particles decouple from the rest of the matter 
at some moment, such that the total density here on 
remains constant. Below the critical temperature T < Tc, 
close to the local minimum the density oscillates around 
the value p = k%K^{T'^ - T'^)/A, being ksnTc ^ 1. The 
function S* in ([B]) has a simple expression, S = sot, with 
So << mc/h in the non-relativistic case. This implies 
that the velocity Vq = 0, also, if there does not exist 
an external force in the system then, = 0. In this 
case the viscosity (dissipative term) of the BEG might in 
fact contain the whole information of the phase transition. 

Finally, to illustrate the previous exposition, we give the 
following example. Suppose that in the system there are 
only condensed and excited particles of the same specie. 
Thus p = Nq/{Ncx + Nf)) (total number of particles), 
where Nq is the number of condensed particles and Nex 
the number of exited particles. Gombining the equations 
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obtained in the theoretical framework lead to the number 
of particles in the condensate in a four dimensional space, 



No 



N 



1.2 ,.27^2 



1 



being N = Nq + Nex , and remembering we are considering 
space and time, note that in this case the exponent 2 
in the critical temperature appears naturally. As always 
this expression shows the dependence of the condensate 
fraction Nq/N as a smooth function of temperature from 
T >Tc down to T = OK. In this case, the finite tempera- 
ture terms are obtained from the one loop corrections of 
the SF density, and are down to be in complete agreement 
with the standard theory, p0l - l22] . 



Observe that only a fraction N / [kl^K^T'^) of the scalar 
particles reaches the ground state at T = 0, this value 
can only be determined experimentally and fits the value 
of the scale k. So in principle we are able to mimic 
the result that in the presence of interactions we have 
A'o < N even at T = 0. The main idea we want to point 
out here is that these phenomena might be equivalent 
for a BEG on earth as for the cosmos, and this might 
follow the previous equations exactly, so this function 
might be tested in the laboratory. As the actually known 
observational evidence favors an open universe in what 
concerns BEG, we therefore need experimental tests to 
measure any new results that can be drown from the 
previous set of thermodynamical equations. If confirmed, 
the phase transition of a BEG can be explained using 
quantum field theory in a straightforward way. 

By adding a J7(l) SB term and temperature contribu- 
tions to the effective Mexican hat potential of a system of 
weakly interacting bosons we have obtained several ther- 
modynamical relations for the nonlinear Schrodinger equa- 
tion ([5]) beginning with KG. We have demonstrated an- 
other way, different from previous works, to obtain the 
equivalence of finite temperature SB and a BEG. 
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